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PREFACE 

This  paper,  Trelimmaiy  Assessment  of  the  Importance  of  Turbulent  Coagulatitxi  in 
the  Kuwaiti  Oil  fnies."  is  part  of  a  larger  Defense  Nuclear  Agency  (DNA)-sponsored 
program  dealing  with  the  Kuwaiti  Oil  Fites  (KC^  of  1991.  This  DNA  program  is  directed 
toward  improving  our  capability  to  predict  the  transport  of  pollutants  arising  from  very 
large  natural  and  manmade  fires. 

The  author  has  received  constructive  comments  and  support  from  numerous 
individuals  in  the  preparation  of  this  paper.  Thanks  are  due  to  Dr.  Ernie  Bauer  (IDA)  for 
pointing  out  the  importance  of  the  fractal  nature  of  aerosols  produced  in  these  fires  and 
the  consequences  of  fiactal  considerations  in  coagulation;  to  Dr.  Ed  Townsley  (IDA), 
Dr.  John  Cockayne  (SAIC),  and  LCDR  Harris  O'Bryant  (DNA)  for  their  critical  review  of 
the  equations  and  mathematical  formalism;  ro  Dr.  Darrel  Baumgardner  (NCAR)  for  helpful 
cormnents  concerning  coagulation  processes  and  the  related  aerosol  size  distribution;  to 
Dr.  Ian  Sykes  (ARAP)  for  pointing  out  the  relaticmship  between  time  averaging  and  plume 
spreading;  to  Drs.  Robert  Oliver  (IDA)  and  Tom  Yonder  Haar  (METSAT)  for  their  general 
review;  and  to  Dr.  Ed  Townsley  and  Mr.  Peter  Kysar  (IDA)  for  their  computational 
support  In  addition,  the  author  would  also  like  to  thank  Mrs.  Sharon  Y.  Wiley  for  her 
outstanding  job  of  technical  typing. 
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EXECUTIVE  SUMMARY 


The  Kuwaiti  Oil  Fires  (KOF)  of  1991  provided  an  opportunity  for  the  Defense 
Nuclear  Agency  (DNA)  to  address  some  key  questims  about  the  atmospheric  transport*  of 
particles  that  are  relevant  to  the  prediction  of  late-time  dust  cloud  nootion.  During  the  fires, 
much  data  were  collected  by  numerous  U.S.  and  foreign  scientific  agencies,  and  now  the 
critical  issue  is  to  assemble  and  organize  a  consistent  set  of  measurements  that  can  be 
evaluated  against  computer  models  to  predia  small  particle  transport.  These  computer 
models  are  required  for  the  prediction  of  long-range  transport  of  nuclear  dust  clouds  and 
for  the  prediction  of  smoke  plumes  from  large  oil,  industrial,  or  urban  Hres  that  might 
affect  electro-t^tical  sensor  performance. 

The  KOF  included  features  that  were  common  with  large  forest  fires,  but  the  KOF 
also  had  other  characteristics  that  made  them  unique.  When  various  long-range  transport 
codes  were  exercised  for  KOF-lUce  problems,  certain  important  deficiencies  were 
encountered.  This  raises  the  following  questions:  Are  the  mathematical  and  physical 
models  of  the  codes  unsatisfactoiy?  Are  they  being  applied  in  situations  for  which  they  are 
not  intended? 

A  key  requirement  for  predicting  particle  transport  is  the  necessity  to  follow  the 
motion  of  some  "tracer”  particles  (Lagrangian  viewpoint).  To  be  good  tracers  for  modeling 
purposes,  the  particles  must  remain  physically  and  chemically  inert  during  the  transport 
process.  This  condition  of  physical  and  chemical  invariance  can  be  termed  a  "frozen"  or 
"aged"  state,  and  it  will  be  reached  when  processes  such  as  condensation  and  evaporation, 
accretion  of  water  vapor  or  other  species,  combustion  and  other  chemical  reactions,  and 
turbulent  coagulation  are  no  longer  active.  In  diis  paper,  we  evaluate  turbulent  coagulation. 
The  other  factors  will  be  addressed  in  subsequent  investigations. 

The  purpose  of  this  paper  is  to  make  a  preliminary  assessment  of  the  conditions 
under  which  turbulent  coagulation  is  or  is  not  a  significant  factor.  Clearly,  if  airborne 
particles  continue  to  coagulate  or  otherwise  grow  or  shrink  in  appreciable  amounts  during 
the  transport  phase,  they  cannot  be  unequivocally  "tagged,"  and  this  potentially  diminishes 


*  In  this  iq)ort,tranqxit  means  both  advection  and  diffusion. 

S-1 


the  validity  of  a  transport  calculadoiL  The  importance  of  this  investigation  is  that  it  defines 
the  nearest  spatial  region  where  turbulent  coagulation  can  be  neglected,  so  that  suitable 
comparisons  can  be  made  between  transport  prediction  models  and  experimental  data 
(notwithstanding  other  aging  processes). 

Specific  reasons  for  evaluating  turbulent  coagulation  in  the  KOF  are  as  follows; 

•  Initial  reported  aerosol  particleAnass  densities  are  extremely  high,  leading  to 
large  coagulation  growth  rates. 

•  In  the  large  Canadian  forest  fire  of  1950,  it  was  found  that  turbulent 
coagulation  was  an  inqxirtant  consideration. 

•  Mass  injection  into  the  atmosphere  is  comparable  to  or  exceeds  that  of  large 
forest  fires. 

•  There  is  a  paucity  of  experimental  data  on  turbulent  coagulation  in  the 
atmosphere,  and  it  may  be  necessary  to  have  a  means  of  interpreting  certain 
KOF  experimental  data  in  this  light 

•  It  is  important  to  understand  the  dependence  of  turbulent  coagulation  on 
windspeed  (as  it  relates  to  tran^xnt  time),  turbulent  kinetic  energy  dissipation 
rate  per  unit  mass,  aerosol  density  and  size  distribution,  and  atmospheric 
diffusion. 

•  The  analytic  representation  of  collection  kernels  for  existing  turbulent 
coagulation  theories  is  based  on  spherical  particles,  but  this  may  not  be  valid 
for  KOF.  Therefore,  it  is  necessary  to  have  a  simple  way  of  estimating  the 
changes  in  coagulation  dynamics  as  a  function  of  aerosol  shape. 

The  approach  is  to  develop  an  analytical  formulation  of  the  spatial  behavior  of 
aerosol  density  as  the  particles  are  canied  along  the  turbulent  wind  field  while 
simultaneously  undergoing  coagulation.  The  KOF  aerosol  problem  is  one  in  which  the 
particle  density  of  the  "supeiplume”  is  formed  by  the  merging  of  hundreds  of  individual 
plumes.  By  example,  Hgure  S-1  shows  three  individual  plumes  merging  into  a  super¬ 
plume. 

In  this  idealized  figure,  we  hypothesize  that  there  will  be  significant  coagulation 
effects  as  the  flames  emerge  from  each  individual  well  since  the  particle  density 
(concentration)  and  mass  density  will  be  extremely  high.  As  the  individual  plumes  spread, 
the  concentration  will  decrease  by  geometric  dilution  caused  atmospheric  turbulence  and 
coagulation.  Ultimately,  we  postulate  that  the  coagulation  process  becomes  insignificant 
when  conq)ated  to  other  processes  due  to  decreasing  concentrations.  The  surfaces  that 


define  this  transition  are  labeled  "Coagulation  Boundary."  The  shaded  areas  behind  these 
boundaries  are  the  regions  where  coagulation  is  assumed  to  be  insignificant 


Figure  S-1.  Geometric  Viewpoint  for  Coaguiation  Considerations 


The  first  step  in  understanding  the  behavior  of  aerosol  concentration  in  the 
superplume  is  to  evaluate  the  contributions  from  an  individual  fire.  Figure  $>2  shows  the 
mathematical  model  for  a  single  plume.  This  model  can  be  applied  to  an  individual  source 
or  to  the  large  plume  provided  that  the  approptiate  initial  conditions  are  present 


AREA  EXPRESSED  AS 
FUNCTION  OF  q 

WIND  SPEED  EXPRESSED 
AS  FUNCTION  OF  q 

AEROSOL  CONCENTRATION 
EXPRESSED  AS  FUNCTION 
OFq 


Figure  8*2.  Modei  tor  individual  Plume 


In  Figure  S-2,  is  the  distance  from  the  source  (oil  well)  where  the  kinetic  energy 
of  the  ejected  oil  is  no  longer  a  factor  in  the  transport  process.  A(qo)  is  the  cross-sectional 
area  associated  with  These  dimensions  arc  estimated  from  visual  observations.  From 
turbulent  diffusion  models,  we  obtain  the  approximation  A(q)  =  (const)  q2  for  q  ^  qo- 
Using  A(q)  in  conjunction  with  mathematical  simplification  of  the  collection  kernels  for 
turbulent  shear  coagulation  and  turbulent  inertial  coagulation  and  an  approximation  for  the 
shape  of  the  aerosol  radius  distribution  based  on  the  large  Canadian  forest  fire  of  1950,  we 
have  been  able  to  derive  an  expression  for  the  aerosol  concentration,  N(q),  as  a  function  of 
the  turbulent  dissipation  rate,  average  windspeed,  and  source  emission  rate.  For  the 
conditions  assumed  in  this  study,  turbulent  inertial  coagulation  is  nearly  two  orders  of 
magnitude  more  efiective  than  turbulent  shear  coagulation. 

The  concentration,  N(q),  along  the  center  line  trajectory  of  Figure  S-2,  is  given  by 

N(q)  =  Nj(q)  rp(q)  ,  (1) 

where 

U(q)A(q)  (2) 

is  the  ordinary  geometric  spreading  and 

rc  =  «p[-i,  (3) 


is  the  part  that  is  due  to  coagulation.  In  these  equations.  So  is  the  particle  emission  rate, 
U(q)  is  the  windspeed  along  q,  A(q)  is  the  cross-sectional  area  of  the  plume,  qo  is  the 
starting  point  for  the  growth  of  the  plume,  and  h  is  a  characteristic  length  associated  with 
turbulent  inertial  coagulation.  It.  is  given  by 


(4) 


where  Rm  is  the  mass  emission  rate,  Pa  is  the  atmospheric  mass  density,  rj  is  a 
characteristic  radius  associated  with  the  aerosol  size  shape  distribution,  Xk  is  the 
Kolmogorov  time  scale  for  turbulence,  Xk  is  the  Kolmogorov  length  scale  for  turbulence, 
and  U  is  the  average  windspeed. 
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Figure  S-3  shows  the  general  behavior  of  Fc  as  a  function  of  q.  The  significant 
features  of  this  curve  are  the  initial  value  of  1.0,  the  asymptotic  behavior  exp(-/i/qo),  and 
the  spatial  variation,  which  behaves  as  exp(-/i/q).  Turbulent  inertial  coagulation  will  not 
have  an  important  effect  on  particle  transport  when  this  factor  does  not  vary  significantly 
over  the  range  of  interest  (i.e.,  the  spatial  region  where  theoretical  models  of  particle 
transput  are  to  be  compared  with  experimental  observations). 


Figure  S*3.  Generic  Behavior  of  Tc  as  a  Function  of  q 


As  shown  in  Figure  S-3,  the  rate  of  spatial  change  of  Fc  decreases  as  the  distance, 
q,  increases.  By  establishing  a  practical  criterion  consistent  with  the  prediction  capabilities 
of  particle  transport  codes,  it  is  possible  to  establish  a  minimum  range,  qmin>  beyond  which 
the  effects  of  coagulation  on  particle  transport  should  not  be  a  concern.  The  criterion  is 


2  0.9  , 


(5) 


and  this  inequality  leads  to  the  minimum  range, 

in  the  domain  of  interest  [where  exp(-/i/qo) «  1].  qmin  is  the  distance  beyond  which  /n  Fc 
changes  by  less  than  1C  percent  over  qmin- 
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Tlte  inapartant  point  is  that  qo  and  the  factors  ccxnprising  h  can  be  determined  from 
renxMe  and  in-situ  measurements  taken  during  the  fires.  Using  estimated  numbers  for  the 
cCTwiitions  of  the  KOF,  we  find  that  /]  can  range  up  to  several  kilometers  and,  in  some 
cases,  can  be  as  large  as  10  kilometers.  On  the  other  hand,  qo  may  only  range  up  to  a  few 
hundred  meters.  At  this  time,  there  remains  a  large  number  of  unknown  factors  in  the 
confutation  of  li  As  the  data  reductitm  and  analysis  from  the  KOF  continue,  improved 
fsriwMitfts  for  li  wiU  become  available,  and  we  will  be  able  to  define  more  precisely  the 
regions  where  a  self-consistent  set  of  measurements  can  be  used  to  validate  particle 
transport  models. 


I.  INTRODUCTION  AND  PURPOSE 


This  paper  is  part  of  an  IDA  Task  about  the  "Characterization  of  Data  and 
Utilization  of  Advanced  Technologies  for  Predicting  the  Atmospheric  Transport  of 
Effluents  frmn  the  Kuwaiti  Oil  Field  Hres,"  being  performed  for  the  Defense  Nuclear 
Agency  (DNA).  The  Kuwaiti  Oil  Hres  (KOF)  of  1991  provided  an  opportunity  for  DNA 
to  address  some  key  questions  about  die  atmospheric  transport  of  particles  that  are  relevant 
to  the  prediction  of  late-time  dust  and  smoke  cloud  motion.  During  the  fires,  much  data 
were  collected  by  numerous  U.S.  and  foreign  scientific  agencies,  and  now  the  critical  issue 
is  to  assemble  and  organize  a  consistent  set  of  measurements  that  can  be  evaluated  against 
computer  codes  to  predict  small  particle  transport  These  computer  codes  are  required  for 
the  prediction  of  long-range  transport  of  nuclear  dust  clouds  and  for  the  prediction  of 
snooke  plumes  fiom  large  oil,  industrial,  or  urban  fires  that  might  affect  electro-optical 
sensor  performance. 

The  KOF  included  features  that  were  common  with  large  forest  fires,  but  the  KOF 
also  had  other  characteristics  that  made  them  unique.  Bauer  (Ref.  1)  has  provided  an  irutial 
assessment  of  the  characterization  of  the  KOF  in  relation  to  various  kinds  of  large  fires  and 
has  identified  certain  key  aspects  that  should  be  studied  further.  However,  the  reviewers  of 
his  document  pointed  out  that  when  various  long-range  tran^>oit  codes  were  exocised  on 
KOF-like  problems,  the  agreement  was  often  not  very  good.  This  raises  the  following 
questions:  Are  the  mathematical  and  physical  rrxxfels  of  die  codes  unsatisfactory?  Are  they 
bang  tqiplied  in  situations  for  which  diey  are  not  intended? 

A  key  requirement  for  predicting  particle  transport  is  the  necessity  to  follow  the 
motion  of  some  "traca"  particles  (Lagrangian  viewpoint).  To  be  good  tracers,  the  particles 
must  remain  physically  and  chemically  inert  during  the  transport  process.  This  condition  of 
physical  and  chemical  invariance  can  be  termed  a  "finozen”  or  "aged"  state  and  aj^lies  to 
those  particles  sufficiently  far  away  fiom  the  sources  of  the  plume.  Some  of  the  processes 
that  should  be  evaluated  to  ensure  that  transport  calculations  for  the  firozen  state  are 
performed  properly  include  condensation  and  evaporation,  accretion  of  water  vapor  or 
other  species,  combustion  and  odier  chemical  reactions.  Brownian  coagulation,  coagulation 
in  laminar  shear  flow,  gravitational  coagulation,  turbulent  shear  coagulation,  and  turbulent 
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inertial  coagulati<xi.  Baumgardner  (Ref.  2)  has  also  suggested  that  turbulent  mixing  may 
have  an  effect  on  aerosol  growdi  under  certain  circumstances.  However,  since  all  of  these 
processes  are  most  active  during  the  initial  stages  of  mass  injection  into  the  atmosidiere  and 
simple  transport  codes  cannot  be  applied  to  a  very  young  plume  in  which  the  tracer  particles 
are  not  yet  "aged,"  it  is  important  to  establish  a  minimum  range  to  which  a  transptnt  code 
can  be  applied  for  particles  whose  mass  no  Itmger  changes  widi  distance  traveled. 

In  this  paper,  we  make  an  initial  assessment  of  the  effects  of  tuii>ulent  shear 
coaguladtHi  and  turbulent  inertial  coagulation  on  aerosol  growth  in  the  KOF.  The  results 
from  this  aiudysis  may  be  useful  in  defining  the  "aged"  regions  that  are  applicable  for 
particle  transport  analysis.  Specific  reasons  for  evaluating  these  turbulent  coagulation 
processes  in  the  KOF  are  as  follows: 

•  Initial  reported  aerosol  paiticleAnass  densities  are  extremely  high,  leading  to 
large  coagulation  growdi  rates. 

•  In  the  Canadian  fire  of  1950,  turbulent  coagulation  was  an  important 
ctmsideradon  (see  Sections  n  and  m). 

•  Mass  ejection  is  comparable  to  or  exceeds  drat  of  large  forest  fires. 

•  There  is  a  paucity  of  experimental  data  on  turbulent  coagulation  in  the 
atmosphere,  and  it  may  be  necessary  to  have  a  means  of  interpreting  certain 
KOF  experimental  data  in  this  light 

•  It  is  important  to  understand  the  dependence  of  turbulent  coagulation  on  wind- 
speed  (as  it  relates  to  transport  time),  turbulent  kinetic  energy  dissipation  rate 
per  unit  mass,  aerosol  density  and  size  distribudon,  and  atmospheric  diffusion. 

•  The  analytic  representation  of  collection  kernels  for  existing  turbulent 
coagulation  theories  is  based  on  spherical  particles.  While  this  may  not  be 
valid  for  KOF,  it  is  necessary  to  have  a  simple  way  of  estimating  the  changes 
in  coagulation  dyruunics  as  afunction  of  aerosol  shape. 

For  noathematical  sin^licity,  nearly  all  analyses  involving  the  transport  of  particles 
from  large  fires  have  assumed  that  these  particles  have  a  spherical  shape.  Since  a  sphere 
has  the  smallest  surface-to-volume  (naass)  ratio,  these  naodels  usually  place  a  lower  bound 
on  the  reaction  rate  per  mass  unit .  However,  unlike  a  forest  fire,  in  which  roughly  half  of 
the  condensed  smoke  particles  are  liquid  and  can  be  described  a  spheres,  an  oil  fire  presents 
a  particular  problem  since  most  of  the  pure  soot  particles  are  Itxig  strands  of  solid  carbon  or 
more  coirq>lex  noolecules  and  cannot  really  be  described  as  spheres. 
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Evans  et  al.  (Ref.  3)  have  examined  the  particle  size  and  shape  distribution  for  oil 
fires  without  interference  fixnn  brine,  water,  nebulized  tar,  or  dusty  ambient  aerosols  and 
have  shown  that  these  smcdce  particles  are  an  agglomeration  of  individual  spherules  that 
exhibit  a  strand-like  behavior.  A  sample  of  their  results  is  shown  in  Figure  M,  which 
indicates  that  large  numbers  of  strand-like  particles  may  have  existed  in  the  many  plumes 
from  the  KOF.  This  particle  configuration  wtxild  have  an  important  effect  on  the  time  scale 
for  kinetic  processes  and  chemical  reactions. 

In  recent  years,  more  attention  has  been  paid  to  the  fractal  description  of  the 
irregular  structures  shown  in  Hgure  I-l  and  on  the  effects  of  such  irregularities  on  free 
molecular  agglomeration  (Ref.  4).  However,  for  the  purposes  of  this  study,  we  will 
assume  a  spherical  shape  for  all  particles  and  sq>ply  a  correction  factor  when  more  accurate 
analysis  becomes  available. 


Figure  1*1.  Electron  Micrograph  of  a  Smoke  Particle  From  a  1.0  m  Diameter 
Murban  Crude  OH  Fire  (Source:  Ref.  3) 


It  should  be  emphasized  that  the  selection  of  the  tuzbulent  coagulation  processes  as 
the  first  of  the  aging  jnocesses  to  be  considered  does  not  imply  that  the  others  have  been 
deemed  to  be  less  importanL  Indeed,  we  plan  to  address  the  other  aging  processes  in  a 
subsequent  publication  in  which  we  will  make  use  of  the  mathematical  formalism 
established  in  diis  study. 

In  Section  II  of  this  paper,  we  evaluate  the  dependence  of  turbulent  shear 
coagulation  and  turbulent  inertial  coagulation  on  aerosol  concentration  and  mass  density, 
size  distribution,  and  turbulence  parameters.  These  considerations  are  then  applied  to  the 
KOF  in  Section  m,  in  which  the  effects  of  windspeed,  atmosphere  turbulence,  and  aerosol 
emission  rate  are  taken  into  account  Concluding  remarks  are  rendered  in  Section  IV. 
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n.  TURBULENT  COAGULATION  PROCESSES 


Because  the  conditions  of  the  KOF  have  not  been  encountered  before,  we  have  a 
limited  data  base  for  predicting  the  spatial  regions  and  conditions  for  which  turbulent 
coagulatitm  may  be  important  It  is  also  unlikely  that  simple  tran^xrrt  codes  will  be  valid 
for  this  type  of  assessment  One  of  the  reasons  for  examining  the  significance  of  turbulent 
coagulatitm  in  die  KOF  is  that  the  initial  reported  particle  concentrations  are  extremely  high. 
Fm  example,  Hobbs  and  Radke  (Ref.  5)  report  particle  concentrations  frequently  exceeding 
lOS  cnr^  within  a  few  kilometers  from  the  fires  and  mass  concentrations  of  the  composite 
plume  ~840  pg  mr^  at  20  km  downwind  for  particles  <3.5  pm  in  diameter.  Johnson  et  aL 
(Ref.  6)  report densities  between  5(X)-1(XX)  pg  m~3  at  distances  greater  than  200  km. 
These  observations  suggest  that  coagulation  processes  could  be  important  over  a  significant 
spatial  regiotL  In  this  sectioi  and  in  Section  m,  we  examine  titis  conjecture. 

Recently,  Pmch  (Ref.  7)  has  examined  the  inq>act  of  turbulent  coagulation  on  the 
size  distribution  of  aerosols  produced  by  large  fires  in  an  attempt  to  understand  the 
phenomenon  of  Blue  Moons.  Porch  used  models  where  initial  mass  concentrations 
for  spherical  smoke  particles  with  radii  less  than  1  pm  ranged  from  5  x  1()~9  to 
5  X  10~S  g  cm~3  (or  equivalently  5()(X)  to  S0,(XX)  pg  m~3).  These  estimates  appear  to 
i^ply  over  the  fire  itself,  represented  in  Porch's  model  by  a  diameter  of  400  km.  In  an 
earlier  study.  Porch,  PCnner,  and  Gillette  (Ref.  8)  conducted  a  numerical  study  of  the 
effects  of  super-pm  particles  <mi  the  coagulation  loss  of  submicron  particles.  The  results  of 
these  studies  have  shown  that  for  particles  with  radii  less  than  0.1  pm  Brownian 
coagulation  is  man  nq)id  than  turbulent  coagulation,  while  the  latter  mechanism  generally 
dominates  for  particle  sizes  of  a  few  microns  (Ref.  9). 

The  basic  model  for  turbulent  coagulation  was  developed  by  Saffriuui  and  Turner 
(Ref.  10),  who  identified  the  two  distinct  coagulation  mechanisms:  turbulent  shear 
coagulation  and  turbulent  inertial  coagulation.  Pruppacher  and  Klett  (Ref.  9)  have 
succinctly  sunomarized  the  salient  features  of  turbulent  coagulation.  The  first  significant 
feature  is  that  turbulent  coagulation  qiplies  to  those  aerosols  of  radius  r  «  where  Xt  is 
the  "Kolmorogov  microscale  length"  (Refs.  9,  10).  As  shown  in  Table  II-I,  this  is 


satisfied  for  the  conditicMis  of  KOF,  in  which  virtually  all  of  the  aerosol  particles  have 
dimenskMis  less  than  die  indicated  values  of 


Table  ll-L  Kolmogorov  MIcrosealaa  as  a  Function  of 
Energy  Dissipation  Rata,  e 


e 

(cm?r^ 

Xk 

(cm) 

tk 

(8) 

Comment 

5 

1.8  xIO-^ 

1.9x10-' 

Appicable  to  straWonn  clouds  where  there  is 
smal  mean  velocity  (Ref.  10) 

100 

8.4  X 10^ 

4.1  X  10-2 

Cbudy  air  (Ref.  9) 

1000 

4.8x10-2 

1.3x10-2 

Relevant  for  condtions  in  turbulent  cumulus 
douds  (Rd.  10) 

4.0  X  10-2 

9.3x10-2 

Relevant  to  early  stages  of  the  plume  in  large 
fires  (Ref.  7) 

2.8x10-2 

4.6  X  10-2 

Exceptional  case  of  turbulent  shear  in  high 
turbulence  (Ref.  8) 

Coagulation  results  from  the  velocity  motion  between  particles.  One  method  of 
achieving  relative  motion  is  through  significant  smaU-scale  velocity  gradients  in  turbulent 
flow.  For  the  length  scale  of  velocity  gradients  that  art  smaller  than  Xk.  the  coagulation 
process  is  similar  to  that  of  die  laminar  shear  flow  case  considered  by  Smoluchowski  and 
described  in  Pruppacher  and  Klen  (Ref.  9).  mi,  m2  are  the  masses  of  the  coagulating 
particles,  and  ri  and  t2  are  their  respective  radii,  and  p  is  the  material  intrinsic  density 
(assumed  die  same  for  bodi/all  particles),  and 


then  die  collection  kernel  for  turbuknt  shear  coagulation  is  given  hy  (Refs.  9, 10): 

“2)  ■  *  '2/  ^  k  •  (2^) 

vdiere  tk  is  the  "Kolmoiogov  time  scale"  (Refs.  9, 10). 

Another  method  for  producing  a  relative  particle  velocity  is  due  to  local  turbulent 
acceleratitHis.  Particles  of  different  mass  respond  differentiy  to  these  accelerations  by 
having  different  viscous  relaxation  times  associated  with  veloci^  equilibration.  This  leads 
to  die  following  expression  for  die  collection  kernel  for  turbulent  inertial  coagulation: 


(2.3) 


K,(n.,.n^)-^^(r,+rj)"  cm’ s'*  . 

In  this  equatitm,  pA  is  the  air  mass  density,  and  and  are  the  Kolnoogorov 
mianscales  for  length  and  time,  which  are  given  by 

/ 

\=(^)  «> 

S  .  (2.5) 

where  e  is  the  dissipation  rate  of  kinetic  enogy  per  unit  mass  (cm2  ^3)  and  v  is  the  kinetic 
viscosity  of  air  (cm2  $-1).  Por  the  altitude  range  used  in  our  study,  we  take  the  U.S. 
Standard  Atmosphere  values  of  pa  -  10~3  g  cnr^  and  v  =  0.172  cn^  s~^ 

As  observed  in  Eqs.  (2.4)  and  (2.S).  the  dissipation  rate  is  a  key  parameter. 
Table  n-I  gives  values  of  Xk  and  Xk  for  selected  values  of  e.  Examination  of  this  table 
shows  that  values  of  e  ranging  from  2000  to  8000  cm2  are  being  considered  in  this 
study. 

Some  of  the  important  mathematical  and  physical  characteristics  of  the  coagulatim 
process  can  be  presented  best  by  considering  tire  time  behavior  of  die  aerosol  distribution  in 
a  ^latially  uniform  mediuoa.  We  let  f(m.  t)  dm  be  the  number  of  particles  per  unit  volume 
lying  in  the  mass  range  between  m  and  m  +  dm.  where  f(m.  t)  is  defined  as  the  mass 
density  function  (MDF).  Using  K(mi.  m2)  as  the  general  description  for  either  of  the 
collection  kernels  [(2.2)  or  (2.3)]  gives  the  following  equation  for  the  time  evolution  of 
f(m.  t): 

2  J  /«(■», *"2 

“1  “2 

-  f(m)  J K(m.  m^^  f(rr^ J  dm^  .  (2.6) 


(2.6) 


The  first  term  in  Eq.  (2.6)  is  the  growth  term.  The  ”1/2”  appearing  in  front  of  the 
integral  is  due  to  the  fact  that  there  is  a  reduction  of  one  particle  when  two  particles 
coagulate.  Inclusion  of  the  Delta  function  S(m  -  mi  -  1D2)  selects  from  all  cnllidnfis  only 
those  in  which  the  noass.  m.  equals  mi  -•>  n^.  It  is  possible  to  integrate  the  growth  term 


over  either  of  the  dummy  variables  mi  or  m2  in  Eq.  (2.6)  and  simplify  the  integral  to  the 
more  fanoiliar  form  (Ref.  7), 

“l“2 

m 

s  ~  jK^m-m^.m^^f^m-m^Jf^m^^dm^  .  (2.7) 

which  is  obtained  by  integration  over  mi.  The  second  term  in  Eq.  (2.6)  is  the  loss  due  to 
coagulation  and  is  die  integradon  over  all  possible  events. 

A  further  condition  that  is  not  explicitly  represented  in  the  expressions  for  Ks  and 
Ki  is  the  existence  of  an  efficiency  factor  associated  with  coagulation.  For  usual  turbulent 
clouds.  Saffinan  and  Turner  (Ref.  10)  suggest  that  the  probability  of  coagulation  may  be 
close  to  unity  when  particles  are  about  the  same  size  but  diminishes  rapidly  for  particles  of 
very  different  size.  The  coagulation  efficiency  of  the  KOF  is  not  clear  since  the  particles 
are  not  water  droplets.  The  issue  of  collection  efficiency  is  difficult  to  address  in  this 
preliminary  assessment,  particularly  since  many  of  the  KOF  particles  may  not  be  sphericaL 
Hence,  for  the  present  we  will  not  include  this  ccmsideiation  in  our  model 

For  this  study,  we  assume  an  efficiency  of  unity  for  all  colliding  particles.  We  also 
assume  that  all  particles  are  spherical  since  the  collection  kernels  have  only  been  modeled 
for  this  case.  A  preliminary  estimate  of  the  collection  efficiency  for  nonspherical  particles 
suggests  that  it  would  be  greater  than  that  for  spherical  particles  of  the  same  mass. 

A  rtxne  precise  evaluation  of  the  collection  kemel(s)  could  perhaps  be  obtained  by 
extending  the  Saffinan  and  Turner  (Ref.  10)  analysis  to  more  realistic  shapes.  However, 
until  such  an  analysis  becomes  available,  we  will  use  the  spherical  noodels.  At  the  end  of 
the  analysis,  it  is  possible  to  apply  a  cotiection  factor  to  account  for  enhancement  in  the 
coagulation  process.  Although  we  do  not  explicitly  include  this  correction  facmr  in  this 
analysis,^  it  is  used  in  the  discussion  of  results. 


1  ActoaDy,  the  combined  effea  of  collisions  and  the  probability  of  particles  stiddng  together  after  the 
collision  are  both  impmtanL  On  one  hand,  the  nse  of  spherical  aerosols  appears  to  underestimate  the 
collection  kernel;  however,  this  is  counterbalanced  in  our  model  by  assuming  that  the  particles  stick 
togedier  in  every  collision.  If  estimates  of  die  relevance  of  turbulent  coagulation  in  KOF  turn  out  to 
be  impocumt  in  long-range  transpat  prediction,  it  miy  be  necessary  to  examine  these  issues  further. 
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Using  Eq.  (2.6),  we  can  make  certain  observadmis  about  the  coagulaticm  dynamics 
to  an  estimate  of  the  relevant  time  constants  as  a  function  of  particle  density  and 
dissipation  rate.  We  first  notB  diat  die  mass  doisity.  M.  wfaeie 


Msjinf(m)dm  , 


(2.8) 


remains  invariant  This  can  be  demcmstrated  by  multiplying  both  sides  of  Eq.  (2.6)  by  m 
and  then  integrating  over  all  mass  space.  Thus,  we  have 


“l“2 

-  J  JmK^ntm^^dm  dm^  . 


“l®2 


JnS(m-m}  - 


(2.9) 


m2)  dm. 


The  first  term  in  Eq.  (2.9)  was  obtained  duough  the  integration 
Since  K(mi,  m2)  and  K(m.  m2)  are  symmetric  kernels,  the  right-hand  side  of  Eq.  (2.9) 
equals  zero  and  dius 

du 

IT-® 

as  expected. 

The  total  number  of  particles  per  cm^  is  defined  as 


(2.10) 


N(t)=Jf(m,t)dm 

and  is  determined  from  Eq.  (2.6)  by  integratitm  over  mass  ^ace.  We  obtain 


(2.11) 


“l“2 


”l"j 

-2  J  • 


(2.12) 


in  in 
1  2 


Improved  insight  into  the  kinetics  of  coagulation  is  gained  by  introducing  a 
normalized  mass  density  functkm  (NMDF),  ^m,  t),  duough  the  equation 


(2.13) 


From  Eq.  (2.1 1),  we  deduce  the  lequiied  nonnaJizaticm  omditicMi 


t)  dm  «  1  . 


(2.14) 


By  introducing  ^m,  t),  we  can  cast  die  analysis  in  terms  of  those  issues  that  are 
attributed  solely  to  the  tmal  aerosol  cmicentration,  N(t),  and  those  issues  that  are  due  to  die 
distributkm  in  mass  as  given  by  ^m.  t).  Substituting  Eq.  (2.13)  into  Eq.  (2.12)  gives 

.  (2.15) 


where 


i1(t)  »  J  J  K^mj.m^J^^mj.t^^^mj.t^dmjdmj  . 


(2.16) 


“i  “2 


vdiere  Nj  is  the  initial  concentration. 

We  will  now  make  an  estimate  of  q(t)  using  the  analytical  forms  of  Ks(mi,  m2) 
and  Kj(mi,  m2)  from  Eqs.  (22)  and  (2.3)  in  omnbination  with  certain  experimental  results 
of  Radke  et  aL  (Ref.  11)  and  as  reported  by  Porch  (Ref.  7),  which  lead  to  an  assessment  of 

^m,t).  Because  Ks  and  Ki  are  easier  to  work  with  in  terms  of  aerosol  radius  instead  of 

mass,  it  is  convenient  in  computing  Ti(t)  to  express  the  mass  distribution  in  terms  of  the 
aerosol  radius  distribution.  Introducing  a  normalized  particle  radius  function  (NPRF), 
Ot  and  ^iplying  the  transformation  law  for  probabili^  frequency  functions  gives 

^m)dm  =  g(r)dr  ,  (2.19) 

where  m  and  r  are  related  by  Eq.  (2.1).  As  seen  in  Eq.  (2.19),  a  knowledge  of  either  or 
g  determines  the  other  differentiating  Eq.  (2.1). 


Fbr  the  turbulent  shear  coagulation  model,  the  expansion  of  (ri  in  Eq.  (2.2) 

gives 

■  ("i  *  *2  *  *2)  •  ® 

Substituting  Eq.  (2.20)  into  Eq.  (2.16)  and  using  the  property  that  Ks  is  symmetric  in  ri 
and  r2  yields 

ilg(t)  »  1.3t"*  (2r^  +  6i^r)  . 

=  2.6t^V(l  +  3r^r/r^)  .  (2.21) 


where 

r^  =  Jr^g(r.t)dr  =  ^  p"*  J  m^m,t)dm  =  ^p"*m  . 
r^  »  J  rg(r,t)dr  , 

?  »  J  rg(r,t)dr  , 

and 

m  s  Jm4Km.t)dm  . 

m  is  the  average  mass  of  all  particles  as  can  be  seen  firom  the  definition 

I  mf(m,t)dm  r 

- - jjjg - . 


(2.22) 

(2.23) 

a24) 

(2.25) 

(2.26) 


Let  us  temporarily  assume  diat  3  r^  r /i^  is  much  less  than  unity  [we  will  partially 
substantiate  this  assutiq)tion  using  Radke’s  data  (Ref.  11)].  Then  1)3(0  is  approximately 
given  by 


2.6  X.  *  ~  p“*  m  . 
k  4jc 


(2.27) 


(2.28) 


Substituting  Eq.  (2.27)  into  Eq.  (2.15)  gives 

dN  0.31 


(Nm)  N 


On  the  other  hand,  since  the  total  tnaj«  density,  M  =  Nm,  remains  constant  during 
the  coagulatitm  process,  turbulent  shear  coagulation  can  be  approximated  by  the  simple 
equation 

aN  _  J1  _ 

at  T  ’  (2.29) 


where 


M 


(2.30) 


The  smallest  values  of  Tg  coirespond  to  small  values  of  tk  (high  levels  of  turbulence)  and 
large  values  of  M.  Using  M  =  5  x  10^  g  cm“3,  p  =  1.0  g  cmr^,  and  tk  =  4.6  x  10-3  s 
(cofresponding  to  e  «  8000  cm^  s-3)  gives  an  e-folding  time  of 


T-  -  2.9x10  s  =  81hr  , 

which  qrpears  quite  long  when  oonqraied  to  the  dme  scales  of  interest  in  the  KOF. 
Let  us  now  return  to  evaluating  the  term 

3i^r/r^  =  A  , 


(2.31) 


(2.32) 


which  first  appeared  in  Eq.  (2.21)  and  wa.<:  subsequently  neglected  in  the  derivation  of 
Eq.  (2.28).  A  may  be  considered  a  "correction  factor"  and  in  the  approximation  of 
Eq.  (2.27)  has  been  neglected.  However,  as  shown  in  Eqs.  (2.22)  to  (2.24),  each  of  the 

entities,  f,  r^,  and  r^ ,  involved  in  die  computation  of  A  is  tune  dependent  since  g(r,  t) 
changes  with  tune.  Essentially,  A(t)  cannot  be  determined  exacdy  without  solving  the 
problem  exacdy.  However,  it  is  interesting  to  explore  the  conditions  under  which  A  is 
neglligible  and  under  which  it  is  greater  than  unity,  the  latter  condition  requiring  modifi- 
catimi  of  Eq.  (2.28). 

An  estimate  of  A  can  be  obtained  using  measurements  of  what  Radke  (Ref.  11) 
calls  die  Number  Size  Distribution.  In  his  notation,  diis  quantity  is  labeled  "dN/dR"  (see 
Figure  II-l)  and  is  d^ned  as  the  number  of  particles  per  unit  radius  per  unit  volume.  Ii 
our  fonnalism,  we  express  diis  quantity  as  N  g(r)  and  have 


(2.33) 


dR 


=  Ng(r)  . 


Particie  racfi  (cm) 

Fiourt  Ihl.  Numbar  Slu  Distribution  of  Particlss  Obssrvttf 
In  a  Forast  Slash  Firs  (Sourca:  Raf.  11) 

The  solid  line  in  Figure  n-1  is  dN/dR  at  the  beginning  phases  of  a  forest  slash  fire. 
The  dashed  line  in  this  figure  is  an  approximation  of  dN/dR,  which  is  necessary  for  die 
approximate  analytical  assessment  used  in  our  study.  The  approximation  of  Figure  n-1 
qiplies  in  die  radius  range  extending  from  a  lower  limit  of  2  x  10-^  cm  to  an  upper  limit  of 
2  X 10-2  cm  and  describes  a  distribution  that  varies  as 

We  should  emphasize  that  die  Number  5ize  Distribution  of  Hgure  II-l  is  talcfn  at  an 
instant  of  time.  However,  during  the  coagulation  process,  g(r,  t)  will  change  with  rime  as 
the  size  distribution  shifts  to  larger  particles  (Ref.  7).  In  addition,  the  g(r)  of  this  figure 
applies  to  a  forest  slash  fire  and  would  likely  be  different  from  the  g(r)  of  die  KOF. 
Nevertheless,  the  use  of  g(r)  from  Hgure  11-1  provides  a  mathematical  basis  for  examining 
the  importance  of  the  radial  shape  function  on  coagulation  dynamics,  and,  therefore,  ura 
assume  that  g(r)  for  the  KCF  is  similar  to  the  result  for  Hgure  II-l. 


The  sensitivity  of  tiie  results  to  the  analytic  representaticm  of  the  aerosol  distribution 
can  be  eiqdained  by  proposing  that  g(r)  for  the  KOF  varies  as  r*  betwe^  n  a  lower  limit, 
and  an  upper  limit,  ixj.  This  form  of  the  apinoximation  for  g(r)  is  chosen  because  of  its 
mathematical  sin^)licity  and  because  it  is  representative  of  aerosol  size  distributions  found 
in  nature  before  the  onset  of  significant  aging  caused  by  coagulation.  For  this  case,  we 
have 

g(r)-Ar«  ,  „  (2.34) 

where  A  is  the  nosmalizatuHi  constant  determined  from  the  requirement  tiiatjg(r)  dr  »  i, 
^tplying  the  nonnalizatkm  condition  gives  ° 

-1 

A  -  (s- 1) rj;*  (l  .  (2.35) 

From  Hgure  II-l,  we  have  it «  2  x  10-6  cm  and  lu  =  2  x  10-2  cm.  The  ratio  TjJtu  equals 
10-4  and  is  tiierefore  negligiUe.  In  this  case,  we  have  the  approximation 

A  =  (s-l)iJ^^  .  (2.36) 


For  expediency,  we  also  assume  that  Eq.  (2.36)  iq)plies  for  the  KOF,  although  this 
assumption  can  readily  be  modified  using  umial  data.  Lastly,  it  should  be  stated  tiiat 
Eq.  (2.34)  implies  a  time-invariant  size  distribution,  which  is  clt.^ly  not  the  case.  Despite 
this  limitation,  Eq.  (2.34)  provides  a  basis  for  understanding  the  kpportant  features  of  the 
size  distributimi  in  regard  to  die  coagulation  process. 


Equation  (2.34)  readily  furnishes  the  eiqnession  for  the  nth  moment: 
?.  Jr»g(r)dr.  . 


(2.37) 


where 

°  ('u'^l)  •  0-38) 

The  moments  r,  r^,  andr^  are  readily  computed  from  Eq.  (2.37).  Inserting  the  resulting 
eiqiressitms  into  Eq.  (2.32)  gives 


n-10 


A 


3(4-s)  (q^-i)(q^-i) 

(2-s)(3.s) 


(2.39) 


We  now  examine  the  dependence  of  A  on  s  in  die  range  where  Q  »  1.  These 
results  are  shown  in  Table  n-2  and  are  cditained  by  taking  the  limit  of  Eq.  (2.39)  for 
Q«r(yiq.»  1.  This  is  the  range  interest  for  forest  fires  but  also  ai^iears  to  be  die  range 
of  interest  fisr  KOF,  based  on  preliminaiy  data  obtained  from  Baumgardner  (Ref.  12).  It 
should  be  noted  that  the  expresskms  in  TaUe  n-2  cannot  be  qiplied  at  the  integer  values  of 
s  «  2, 3, 4.  In  order  to  cmiqiute  A  at  diese  values,  the  limiting  form  of  Eq.  (2.39)  at  s  «  2, 
3, 4,  respecdvely,  must  be  taken  using  UHopital's  rule.  Hgure  n-2  provides  a  broader 
perqiecdve  of  die  dependence  of  A  (m  s  and  Q. 


Table  11-2.  Dependence  of  A  on  a 


8 

A 

1  <8<2 

3(4-81 
(2-8)  (3-8) 

r 

2<8<3 

3(4:8)  / 

(8-2)  (3H1)  \ 

m 

3(4-8)  / 

(8-2)  (8-3)  \ 

4-8 

4<8 

3(8-4) 

(8-2)  (8-3) 

Examination  of  Table  n-2  shows  that  the  conection  factor  A  is  negligible  for  s  >  2 
since  «  1.  On  the  other  hand,  in  the  range  1  <  s  <  2,  the  variation  of  A  with  r^At 
shows  that  the  correction  factor  can  be  several  orders  of  magnitiMlf*  greater  than  unity.  This 
is  also  indicated  in  Figure  n-2.  Should  g(r)  vary  as  r~B  with  1  <  s  <  2,  the  time  constant, 
Ts,  in  Eq.  (2.30)  could  be  reduced  by  several  orders  of  magnitude,  rendering  turbulent 
shear  coagulation  a  very  important  consideration  since  the  effective  time  scale  would  be 
reduced  to  ^iproximately  1  hour.  It  is  inqiortant  to  note  that  significant  coagulation  effects 
involving  time  scales  under  1  hour  and  for  particles  greater  than  1  pm  in  radius  have  been 


conqnited  for  die  large  forest  fire  of  1950  using  Radke's  results  in  Hgure  II-l 

(Ref.  7). 


Figure  ll^a.  A  es  a  Function  of  a  With  Q  as  a  Parameter 


Up  to  this  point,  we  have  cmly  considered  the  ocmtribution  firom  turbulent  shear 
coagulaticm  through  use  of  the  oollectitMi  kernel  given  1^  Eq.  (2.2).  We  now  consider  die 
contribution  from  turbulent  inertial  coagulation  in  Eq.  (2.3).  Using  Eq.  (2.16),  the 
coneqmnding  parameter,  Tii,  is  given  by 

r,  u 


n^iere 


2k 

9 


B 


(2.41) 


n-12 


In  contrast  to  the  case  for  turbulent  shear  coagulatitm,  the  double  integratitm  over 
ri,  12  is  not  sqMiable  due  to  the  nonanalydc  function  ~  •  However,  it  is  possible  to 

obtain  a  lower  bound  on  in  dw  following  way.  Ccm^der  the  integrals  Ii  and  I2,  where 

'1  '2 
and 

*2  -  J  J  ('■  *  '2)"  ('1  -'2/  *(^)  *(y  *1  *2  • 

'1  '2 

1]  ^>peais  in  Eq.  (2.40).  I2  is  similar  and  has  the  essential  propeity  required  of  I  j,  namely, 
ensuring  diat  no  contribution  conies  from  the  domain  where  rj »  r2. 

We  now  show  dial  I|  >  I2.  Let 

a-  lt^-i^|.(t,+rj)lr,-rjl  ^44, 

and 

'’“('i-'2)*“('.-'2)('r'2)  e-'**) 

and  ocmsider  die  difference 
a  -a*-b^ 

’('i*'2f("i-'2f-Ci-y*  ■('r'2f  ^  °  •  e--*® 

Sincea,bare>0,anda22b?,  wectmcludethata^b.  Therefore, 

•  (2-47) 

A  lower  bound  on  qi,  defined  as  i)*,  is  thus  given  by 

Tf  «  B  .  (2.48) 


Use  of  provides  a  minimum  rate  of  coagulation  by  inertial  turbulence.  The  integration 


of  I2  yields 


n-13 


Hientiod^CT^)  is  deduced  from  Eq.  (237)  to  yield 
„  (?f  s-s  (q"-0' 

/  (3~s)^  (q^-i) 


(2.49) 


(2.50) 


Table  n-3  shows  die  dependence  of  6  <xi  s  for  the  range  of  interest  ii^Al  «  Q  »  1.  and 
Figure  n-3  shows  6  as  a  function  of  s  with  Q  as  a  parameter.  It  should  be  noted  that  the 
expressitxis  in  Table  n-3  cannot  be  applied  at  the  integer  values  of  s  »  3,  s.  At  these 
values  of  Se  0  must  be  cooqnited  fimn  Eq.  (2.50)  using  Lllopital's  rule. 


Table  11-3.  Oapandance  of  0  on  a 


a 

_ 8 _ 

1  <8<3 

{3-S)*  * 

•-1 

3<8<5 

(8-3) 

r 

5<8 

8-5 

(8-3)* 

Since  the  anticipated  range  of  interest  is  firtmi  s  «  3  to  s  >  5  and  for  Q  »  1,  we 
conclude  from  both  Table  n-3  and  Hgure  n-3  that  0  can  be  neglected  in  conq;>arison  to 
unity  and  thus  use  die  appr>;xkr>ation 


2r^ 


(5-s) 


(2.51) 
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Figur*  11-3.  6  M  a  Function  of  a  With  Q  aa  a  Parametar 


Using  the  results  fircmi  Table  n-2  and  neglecting  the  contribution  firom  A  for  s  >  2 
provide  a  comparison  of  the  relative  ctmtributitHi  for  shear  and  inertial  coagulation.  For 
s  >  2,  we  have 


and 


(2.52) 


is.  -E_ 

»  Pa 


» 


(2.53) 


and  die  ratio  is  given  by 


(q^-0 


(2.54) 


In  the  subsequent  calculations,  we  use  p  =  1  g  cmr^  for  the  aerosol  particles,  and 
Pa  *  10-3  g  ciir3.  There  is  drainadc  variation  of  tiVt]  as  a  function  of  s.  This  quantity 

1  w 


n-15 


can  vary  over  several  orders  of  magnitude,  depending  on  the  ratio  ixj/ix.  (as  long  as 
r^yitL  »  1).  For  exanq)le,  when  2  <  s  <  4,  Eq.  (2^)  reduces  to 


(540) 


(4-s)  ^ 
(5-s)  3^  ’ 


(2.55) 


while  for  s  >  5  we  have 


(2.56) 


The  ratio  between  Eqs.  (2.55)  and  (2.56)  is  approximately  given  by  ruAL>  which  may 
range  several  orders  of  magnitude  (e.g.,  ru^  « 10^  in  Hgure  n-1). 


It  is  interesting  to  examine  in  dre  range  2  <  s  <  4  for  a  representative  large 

value  of  disidpatifin  rate.  Using  the  value  ^  «  4  x  10~2  cm  from  Table  n>I,  corresponding 
toe«2000,anditj«2x  10-2anddievaltress3deducedfromHguren-l  gives 


i  ^ 


135  . 


(2.57) 


If  were  to  be  used  in  the  formalism  of  Eqs.  (2.28)  to  (2.30),  the  time  constant  would  be 

Ti  instead  of  Ts  and  would  be  givoi  by 

T*  2  o  X  10^  3 

Tj  -  =  135  *  2.1  X  10  s  *  35  min  .  (2.58) 

The  results  show  diat  inertial  coagulation  is  the  dominant  mechanism  for  particle 
frequracy  functions  that  ccmtain  a  substantial  fraction  of  particles  greater  than  1  pm.  This 
conclusitm  is  also  cmisistent  with  die  general  observatioos  made  by  Pruppacher  and  Klett 
(Ref.  9). 


in.  APPLICATION  TO  KOF 


The  purpose  of  this  section  is  to  examine  the  geometric  boundary(ies)  of  KOF 
where  turbulent  coagulation  no  Itmger  contributes  to  the  aging  of  aerosols.  When 
combined  with  the  results  of  other  aerosol  mass  changing  processes,  this  analysis  will 
contribute  to  defining  the  overall  spatial  region(s)  of  KOF  where  aerosols  no  longer 
experience  mass  changing.  These  are  the  only  regimis  where  it  is  meaningful  to  compare 
the  concentration  prediction  ci^Mbility  of  cuirent-genexation  particle  tran^x»t  codes  (which 
do  not  include  aging  effects)  to  experimental  results. 

Figure  in*l  depicts  the  KOF  plume  behavior  for  turbulent  coagulation  consid¬ 
erations.  This  figure  portrays  three  individual  plumes  merging  into  a  superplume.  The 
actual  KOF  aerosol  problem  is  one  in  which  the  particle  concentration  of  a  superplume  is 
forged  by  the  merging  of  hundreds  of  individual  plumes. 


Figure  III-1.  Geometric  Viewpoint  for  Coagulation  Considerations 


In  this  figure,  we  hypothesize  tiiat  there  will  be  significant  coagulation  effects  as  the 
flames  emerge  from  each  individual  well  since  the  concentration/mass  density  will  be 
extronely  high.  As  the  individual  plumes  spread,  the  concentrations  will  decrease  by 
geometric  dilution  caused  by  atnoospheric  turbulence  and  coagulation.  Ultimately,  we 
propose  that  the  coagulation  process  beconxs  insignificant  due  to  decreasing  particle 

ra-1 


concentrations.  The  surfaces  that  define  this  transition  are  labeled  "Coagulation 
Boundary."  The  shaded  areas  beyond  these  boundaries  are  the  regions  whoe  ctuigulatitHi 
is  to  be  insignificanL 

Hgure  m-l  suggests  the  possibility  that  the  mo’ging  plumes  may  again  produce 
conditions  for  coagulation;  hence,  anotho’  Coagulation  Boundary  is  shown  in  the 
superplume.  However.  Sykes  (Ref.  13)  has  questioned  the  possibility  of  a  seccmdary 
coagulation  boundary  for  the  merged  plumes  on  the  grouiuis  that  this  could  not  occur 
without  some  mechanism  for  ctmcentrating  particles  locally.  This  issue  will  be  resolved 
when  detailed  and  reliable  data  are  availatde  fixxn  the  field  ejqKtimmt 

The  specific  objective  of  this  section  is  to  develop  a  mathematical  model  for 
predicting  the  downwind  aerosol  density  for  eitiier  an  individual  fire  or  a  noerged  plume 
and,  firom  diis  nxxlel,  to  determine  the  Coagulaticm  Boundaries  shown  in  Hgure  IQ-l.  The 
aiudytical  model  is  devel<^)ed  1^  axnbining  the  turbulent  coagulation  kinetics  described  in 
Section  n  with  particle  transport  and  incoipcnating  advection  and  turbulent  difiusioiL 

In  Section  n,  we  showed  that  the  aerosol  concentration  for  a  spatially  uniform 
medium  is  given  by 

^  =  ,  (3.1) 

where  11(0  can  be  written  as 
Ti(t)  =  J  J 

^  '2 

Equation  (3.2)  is  the  same  as  Eq.  (2.16)  but  is  now  written  in  terms  of  the  radius  variables. 
Summarizing  the  results  from  Section  n,  we  have  the  following  for  turbulent  shear 
coagulation: 

*  2.6  x;;*?  (1+A)  .  (3.3) 

where 

A  =  3r2  r/r^  .  (3.4) 

For  turbulent  inertial  coagulation^  we  have 


where 


e  > 


(3.6) 


It  should  be  noted  that  Eq.  (3.3)  is  a  restatement  of  Eq.  (2.21)  and  Eq.  (3.5)  is  a 
restatement  of  Eq.  (2.48)  using  the  exact  expression  of  I2  of  Eq.  (2.49)  instead  of 
neglecting  6. 

As  imlicated  in  Eqs.  (3.3)  to  (3.6).  a  knowledge  of  the  first  four  moments  of 

g(r,  t).  r,  and  r^,  precisely  determines  the  rime  evolution  of  the  aerosol  density,  b 
should  be  stressed  that  riiese  quantities  iqiply  to  any  fiincri(»al  foim  of  g(r,  t).  Thus,  any 
measurements  of  aerosol  properties  that  could  furnish  information  regarding  die  Erst  four 
moments  would  be  sufficient  to  determine  the  aerosol  dynamics.  The  situation  becomes 
sinqrlified,  scnnewhat,  when  we  have  further  knowledge  of  A  and  6. 


Although  it  is  not  explicitly  stated,  Eq.  (3.1)  applies  to  the  case  in  which  an  initial 
number  of  particles.  No.  are  insetted  into  the  system  at  rime  t »  0.  The  concentration  at 
subsequent  times  is  given  by 


N 


N, 


1+NqT/2 


(3.7) 


where 

t 

t  =  J  Ti(t')dt'  .  (3.8) 

o 

Equations  (3.7)  and  (3.8)  are  deceptively  sinqrle  because  they  assume  a  knowledge 
of  g(r,  t).  As  we  mentioned  at  the  outset  of  Sectimi  n,  g(r,  t)  cannot  be  determined  without 
solving  the  entire  problent  However,  estimates  of  the  significant  times  scales  involved  in 
the  coagulation  process  can  be  obtained  using  approximate  forms  of  g(r,  t)  or  real-time 
measurements  of  g(r,  t).  It  is  through  the  latter  that  we  consider  turbulent  coagulation  in 
die  context  of  die  KOF. 

In  the  remainder  of  this  section,  we  will  calculate  the  aerosol  concentration  as  the 
particles  emerge  from  a  single  well.  The  mathematical  model  developed  for  this  case  is 
direcdy  applicable  to  the  superplume,  with  die  appmpnatt  change  of  coordinates. 


in-3 


When  the  small-scale  tmbulent  fluctuations  are  averaged  out,  the  Euleiian  equation 
for  the  aerosols  is  given  by  the  extension  of  Eq.  (2.6)  to  now  include  convection.  We  have 

-f(m)  jK(ni.m2)f(m2)dm2  ,  (3.9) 

“*2 


where  U  is  the  average  velocity  and  V  is  the  gradient  operator.  The  right-hand  side  of 
Eq.  (3.9)  is  the  same  as  that  of  Eq.  (2.6).  For  brevity,  we  have  suppressed  the  explicit 
dependence  of  f  on  the  space  variable,  1^  and  on  time,  t.  For  example,  f(m)  is  really 
f(m.  K,  t). 


For  an  individual  plume,  Eq.  (3.9)  is  solved  by  matching  the  solution  to  the  source 
emission  at  the  origin  Oocation  of  the  burning  well).  If  the  source  emission  does  not  vary 
with  time,  Eq.  (3.9)  will  also  be  time  invariant  and,  hence. 


df(m,  ^  t) 
dt 


(3.10) 


Ibe  space  dependence  of  f(m),  ^  is  found  by  constructing  a  simplified  geometric 
model  of  the  plunoe  as  shown  in  Hgure  in-2.  In  this  figure,  we  approximate  the  plume  as 
a  tube  of  variable  cross  section,  A(q),  in  iriiich  all  the  particles  move  with  average  velociQr, 
U,  normal  to  the  surface,  A.  The  concentration  is  assumed  to  be  uniform  across  die  area. 
The  coOTdinate  system  is  constructed  so  that  the  q-axis  is  always  directed  along  the 
direction  of  U.  If  |  IT|  is  the  magnitude  of  U  along  the  q-axis,  the  divergence  term  in 
Eq.  (3.9)  becomes 

V  .  [P«ni.  E)]  -  ^  [U(q)  Hin,  q)]  (3.11) 


since,  by  construction,  f(m,  q)  does  not  vary  in  directions  normal  to  q. 

We  now  let 

f(m,q)  =  N(q)9(ni.q)  ,  (3.12) 

where  N(q)  is  the  aerosol  concentration  (cmr^)  and  ^(m,  q)  is  the  NMDF  (expressed  as  a 
function  of  the  space  coordinate,  q,  instead  of  a  function  of  time). 
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A(q)  «  AREA  EXPRESSED  AS 
FUNCTION  OFq 

U(q)  «  WmO  SPEED  EXPRESSED 
AS  FUNCTION  OF  q 

N(q)  >  AEROSOL  CONCENTRATION 
EXPRESSED  AS  FUNCTION 
OFq 


Figure  IlhS.  Modal  for  individual  Pluma 


Using  Eqs.  (3.10)  and  (3.11)  in  Eq.  (3.9)  and  then  integrating  over  mass  space 
using  the  techniques  developed  in  Section  n  g^ves 

(ANU)  ■  .  (3.13) 

where  the  q-dqiendent  time  "constanL"  T(q),  is  given  by 

“  N(q)^Tl(q) 
and 

il(q) .  J  J  K[m,(r,)  g(t,,q)  g(rj,q)  dr^dr^  .  (3.15) 

'2 

The  expressions  for  ti  corresponding  to  turbuloit  shear  and  turbulent  inertial  coagulation 
are  given  Eqs.  (3.3)  and  (3.5),  respectively.  As  in  Section  n,  we  have  also  introduced 
an  NPRF,  g(r,  q),  through  the  transformation 

^m,  q)  dq  »  g(r,  q)  dr  .  (3.16) 

Further  simplificatitxs  is  possible  using  the  conservation  of  mass.  Multiplying 
Eq.  (3.9)  1^  m  and  then  integrating  over  mass  space  and  area  gives 


^(U(q)M(<i)A(q))  -  0  .  (3.17) 

where  M(q)  is  die  mass  per  cm^.  Equadoi  (3.17)  is  analogous  to  Eq.  (2.10).  is  the 
total  mass  emission  rate  (g  s~*)  from  a  burning  well,  and  as  long  as  the  flux  tube  of 
Figure  in-2  encompasses  all  the  mass,  we  must  have  the  equality 

U(q)M(q)A(q)  -  R„  .  (3.18) 

where  Rm  is  the  mass  emission  rate  from  the  well  expressed  in  g  s~l. 

The  next  step  in  determining  N(q)  is  to  evaluate  A(q).  We  assuma  a  circular  plume 
with  radius,  Ro.  and  area, 

A(q)  »  xRj(q)  .  (3.19) 


Rq  is  qiproximated  by  die  equatkm 

t 

Ro  *  2  JiXtOdt'  ,  (3.20) 

0 

where  D  is  the  turbulent  diffusion  ccmstant  evaluated  along  the  path.  Following  basic 
thetnedcal  arguments  discussed  by  Tennekes  and  Lumley  (Ref.  14),  we  assume  the 
relationship 

D(t)-R^(t)U(t)  ,  (3.21) 

which  applies  for  times  Itmg  enough  for  many  independent  eddies  to  establish  fully 
developed  tutbulence. 

^th  diis  assuiiq>ti(Mi,Eq.  (3.20)  can  be  written  as  a  precise  equation: 

DcaRU  .  (3.22) 

where  "a"  is  an  experimentally  deduced  constant,  as  yet  undetermined.  Inserting 
Eq.  (3.21)  into  Eq.  (3.20)  and  differentiating  the  latter  gives 

dR- 

*  aU(t)  .  (3.23) 


Integrating  Eq.  (3.23)  along  die  path  yields 


(3.24) 


t 

Ro^q)  »  a  J  U(t)dt  «  aq  . 

o 

For  the  turbulent  diffusion  case  considered,  the  radius  of  the  plume  is  observed  to  be 
proportional  to  the  distance  traveled.^ 

In  the  regions  described  by  Eq.  (3.21),  it  is  recognized  that  Eq.  (3.24)  is  an 
approximation  since  the  diffiisicxi  rates  far  hmizmital  and  vertical  diffusion  are  different 
For  tl»  KOF,  it  has  been  observed  that  the  plumes  are  typically  constrained  between  1  and 
3  km  in  diickness.  When  the  top  of  plume  slows  its  rise  and  stabilizes  from  the  initial 
buoyancy  effect  the  padi  variable  q  can  be  replaced  by  die  horizontal  axis  in  the  direction 
of  flow,  x;  the  vertical  dimensicm  can  be  replaced  by  a  constant  value.  Ah  -2  km;  and  the 
horizontal  dimension  normal  to  the  flow  can  be  allowed  to  grow  at  a  turbulent  diffiisicm  rate 
given  by 

y  =  bx.  (3.25) 

where  b  is  another  experimentally  deduced  quantity.  If  it  is  desired,  die  approximadon  for 
A(q)  could  be  improved  by  using  real-time  KOF  eiqxnimental  data. 

However,  until  such  KOF  experimental  data  become  available,  we  will  use  the 
circular  model  described  by  Eq.  (324)  with  the  constant  "a"  determined  from  previous 
experimental  data.  Figure  III'3  shows  the  horizmitai  mean  cloud  half-width,  Oy(x),  for  a 
stack  plume  as  a  function  of  downwind  distance,  x,  from  the  source  OE^efs.  IS,  16). 
Lines  A  through  F  correspond  to  different  levels  of  stability.  Taking  the  "average"  of  the 
curves  gives  the  approximate  relationship 

‘0.1x  .  (3.26) 

Identifying  the  Gy  and  x  of  Eq.  (3.26)  with  the  Ro  and  q.  respectively,  of 


Eq.  (3.24)  gives 

R(q)»0,lq 

(3.27) 

and 

A(q)  =  JiR^  s  3.14  X 10  ^  q^  cm^ 

(3.28) 

^  L  Sykes  (Ref.  13)  has  pdnted  oat  that  die  linear  behavior  of  Eq.  (3.24)  is  an  "early  time”  result  but 
may  well  extend  to  dimensioiis  near  global  scales.  For  very  lo^  transport  times,  die  horizontal 
growth  will  vary  with  the  square  root  of  the  distance. 
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At  this  point,  it  is  interesting  to  estimate  the  mass  densities  calculated  from 
Eq.  (3.30)  for  an  "average"  burning  well  and  compare  them  with  the  values  used  by  Porch 


(Ref.  7).  If  Rm  is  the  aerosol  emisaou  rate  for  all  the  burning  wells  and  W  is  the  total 
number  of  burning  wells,  then  die  average  mass  emission  rate  per  well,  Rmt  is 

^  ^  ■  (3*31) 


We  do  not  currently  have  information  regarding  the  maximum  emissicm  rates  for  the 
largest  individual  fires,  and  there  appears  to  be  some  uncertainty  regarding  both  the  total 
emission  rate,  Rm>  and  the  average  emission  rate.  Bauer  (Ref.  1)  has  estimated  Rm  to  be 
4S,(XX)  metric  tons  per  day,  while  Hobbs  and  Radke  (Ref.  S)  have  estimated  Rm  to  be 
34(X)  metric  tons  per  day.  However,  it  is  not  clear  whether  diese  numbers  apply  to  the 
same  number  of  burning  wells  and  other  similar  conditions.  For  lack  of  more  con^lete 
information,  we  take  the  average  of  the  two  references  and  use  an  emission  rate  of  24,000 
metric  tons  per  day,  which  converts  to  Rm  »  2.8  x  10^  g  s~^  If  we  assume  that  3(X)  wells 
are  buming,^  the  average  emission  rate  per  well  is 


^_2Wl£l..  930gs-> 


(3.32) 


We  may  anticipate  that  all  wells  will  not  be  burning  at  the  same  rate  and  that  some  wells 
may  have  an  emission  rate  far  exceeding  that  given  by  Eq.  (3.32). 

In  order  to  examine  a  hypothetical  "worst-case”  scenario,  we  use  a  value  of 
Rg, «  9.3  X  103  g  s-1  in  the  calculations,  which  is  an  order  of  magnitude  higher  than  the 
average  of  Eq.  (3.32).  At  a  distance  of  1  kilometer  from  the  well  (qjc  =  1),  with  a 
windspeed  of  10  km/hr  (Uw  =  10)  and  well  emission  rate  Rm  s  R^  »  9.3  x  103  g  s~l,  the 
mass  density,  M,  equals  1.06  x  10-7  g  cnr3,  which  is  a  factor  of  2.12  higher  than  the 
maximum  value  of  S  x  10-^  g  cnr^  used  by  Porch  (Ref.  7).  For  die  same  values  of  q^  and 
Uw,  but  with  Rm  *  R^  =  930  g  s-1,  M  =  1.06  x  10-8  g  cmr3,  which  is  consistent  with 
Porch's  (Ref.  7)  range  of  interest,  5  x  10-9  to  5  x  10-8  g  cm-3.  Using  Uw  =  10, 
R^j  =  R^,,  but  this  time  using  q^  =  10,  gives  M  =  1.06  x  10-9  g  cmr3,  which  is  a  factor  of 
4.7  lower  than  Porch's  (Ref.  7)  minimum  of  5  x  10-9  g  cm-3.  The  net  result  of  this 
discussion  is  the  reassuring  observation  that  the  values  of  M  calculated  for  a  realistic  range 
of  conditions  are  consistent  with  the  values  previously  considered  in  a  large  forest  fire. 

We  now  return  to  the  main  thrust  of  this  section,  which  is  the  determination  of 
aerosol  density,  N(q),  from  Eq.  (3.13).  Using  the  result 


7  In  ifaeMay/liine  1991  tune  period,  it  is  estimated  that  between  SOO-600  wells  woe  burning. 
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U  »  N(q)m(q)  -  ^  N(q)  ?(i) 


in  Eq.  (3.18)  gives 


4jtpr^AU 


(3.33) 


(3.34) 


Substituting  Eq.  (3.34),  in  conjuncti(»  with  the  expressions  for  T](q)  obtained  from 
Eqs.  (3.3)  and  (3.5),  respectively,  in  Eq.  (3.14)  gives  the  following  expressions  for  T(q) 
fac  turbuUnt  shear  (Ts)  coagulaticm  and  furbidenr  inertia/ (7/)  ooagulaticm: 

8jtp  AUt. 

V‘»"'75r^  •  (3-35) 


where 


fs=l+A  . 


Tj(q) 


where 


r^(i-e) 


(3.36) 


(3.37) 


(3.38) 


Using  either  of  the  expressions  Ts(q)  or  Ti(q)  for  T(q)  in  Eq.  (3.13)  yields  the 
doisity  along  the  path.  Solving  Eq.  (3.13)  for  "ANU”  gives 


/A„N  U„\ 

r,  . 


where 


(3.39) 


(3.40) 


U(qOT(qO  • 


(3.41) 
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The  initial  values  Ao.  No.  and  Uo  are  taken  at  the  location  qo.  where  the 
approximaticm  of  Eq.  (327)  is  valid  or,  equivalently,  where  the  initial  jetting  action  and 
buoyancy  effects  of  the  well  are  dissipated.  That  is.  our  turlnilent  diffusion  nx>del  cmly 
ctmsiders  those  effects  attributed  to  usual  atmospheric  processes. 

For  the  purposes  of  this  analysis,  we  must  make  an  assumption  about  the  way 
IKq*).  and  TCq')  are  treated  in  the  integration  of  Eq.  (3.41).  Until  experimental  data 
become  available,  we  ^>proximate  these  quantities  as  "path-averaged"  constants.  On  the 
other  we  do  include  the  variation  of  A(q)  from  Eq.  (327).  The  results  for  turbulent 
shear  (S)  coagulation  and  turbulent  inertial  (I)  coagulation  are  as  follows: 


where 


*  ■  u"t.  p 


(3.43) 


and 


(3.44) 


where 


(3.45) 


In  die  absence  of  ooagoladcm,  die  spatial  distribution  of  density  is  given  by 

A„N„U^ 

*  U(q)A(q)  •  (3-46) 

This  spreading  behavior  of  the  plume  is  attributed  to  atmospheric  diffusion.  Ni(q)  can  be 
estimated  from  the  mass  ejection  using  die  equation 

MamN  (3.47) 

in  Eq.  (3.18).  The  result  is 

N,(q)  a  rL--.  ”  . 

mU(q)A(q) 


(3.48) 


Ni(q)  is  die  particle  density  that  would  exist  in  the  absence  of  coagulation.  Since  Pc  ^  1, 
coagulatkn  reduces  the  aerosol  concentration  beyond  ordinary  geomettic  dilution.  Coagu¬ 
lation  is  quantitatively  manifested  by  die  parameters  1$  or  /j,  which  pertain  to  turbulent 
shear  coagulation  and  turbulent  inertial  coagulation,  respectively.  If  we  tenqxnarily  let 
/ « /$  or  /i,  dien  the  coagulation  contribution  to  density  can  be  written  as 

which  is  observed  to  be  less  than  or  equal  to  unity  for  the  prescribed  range  q  ^  qo-  The 
density  at  qo  is 

^(^)*^l(qo)“  SU(q^)A(^  •  (3.50) 

It  is  difficult  to  theoretically  predict  qo  since  the  processes  describing  the  aerosol 
behavior  as  die  oil  emerges  from  the  well  are  extremely  complex.  However,  estimates  of 
qo  can  be  obtained  from  visual  observations.  While  a  knowledge  of  qo  is  required  for 
determining  the  absolute  level  of  aerosol  ctmcentration  for  q  ^  qo,  it  does  not  rqipear  to  be  a 
critical  factor  in  determining  whether  to  consider  coagulation  effects  when  interpreting 
transport  behavior.  What  is  important  is  the  spatial  variation  of  Pc  as  a  function  of  q  when 
conqiared  to  other  transport  processes. 

It  is  now  worthwhile  to  examine  certain  features  of  N(q),  including  predictions  of 
aerosol  density  close  to  the  source  and  the  apparent  limiting  effects  of  coagulation.  Based 
on  crude  visual  interpretations  of  the  KOF  photographs,  we  estimate  the  initiation  of 
"ctmventional"  atmospheric  diffusion  processes  to  begin  between  0.1  and  1.0  kilometers 
from  an  individual  burning  well.  For  orientation,  we  compute  Ni(qo)  at  qo  =  10^  cm 

(1  km)  for  a  windspeed  Uw  =  10  km  hr~i,  for  Rm  =  R^,  and  for  the  average  mass, 

m  ^m  «  (41CP  r^/3}  «  6.7  x  10  g  j,  as  determined  from  the  size  distribution  defined 

Eq.  (A.6)  of  Appendix  A.  Inserting  these  values  in  Eq.  (3.48)  gives  Ni  (at  1  km) » 
1.6  X 10^  particles  per  cm^.  This  value  of  density  appears  to  be  consistent  with  the  results 
of  Hobbs  and  Radke  (Ref.  5)  at  close-in  ranges. 

As  a  matter  of  theoretical  interest,  it  is  worthwhile  to  examine  the  behavior  of 
aerosol  density  as  a  function  of  emission  rate,  Rn.  By  the  way  of  example,  we  evaluate 
the  asymptotic  behavior  of  N(q).  In  the  limit  where  q  » /  and  q  ^  qo,  the  coagulation 
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functim  approaches  the  constant  •  exp  (-//Qo)*  Eqs.  (3.43)  and  (3.45),  we 
note  that  bodi/s  and /i  are  proporticHud  to  Rm.  Therefore, 

/-aRo,.  (3  51) 

where  a  is  the  prcqxirticMuility  con^ant  [for  example,  a  =  (5.3  in  case  of 

turbulent  inertial  coagulatitw].  Substituting  Eq.  (3.51)  into  yields  the  following  result 
for  die  asynqitotic  behavior 

°  (aUWAw)  “K-V*-]  •  (3.52) 

where 

X  =  qya  .  (3.53) 

Keeping  all  other  factors  the  same,  Nu  has  a  maximum  with  respect  to  Rm.  This 
ma«ifnum  results  from  the  competing  factors  of  increased  emission  and  reduction  by 
coagulation.  Nas  maximizes  at  Rm  1,  with  the  maximum  value  given  by 

,  ,  (0.37)  X 

{"tsjau  mU(q)A(q)  ' 

As  we  indicated  earlier,  the  main  issue  regarding  the  relevance  of  turbulent 
coagulaticm  is  to  exaixtine  its  effect  on  the  interpretation  of  transport  behavior.  More 
specifically,  we  want  to  ensure  that  N(q)  does  not  change  appreciably  due  to  turbulent 
coagulation  over  distances  diat  are  significant  in  the  interpretation  of  transport  This  means 
that  we  would  like  the  function  Pc  to  remain  essentially  constant  This  condition  is 
examined  most  easily  by  plotting  Pc  as  a  function  of  q  with  qo  and  /  as  parameters. 

To  obtain  an  estimate  of  values  of  /  that  are  of  interest  we  calculate  Is  and  /j  for  die 
following  set  of  parameters^: 


ID 

*R^=  9.3x  longs’* 

(a) 

P 

s  1  gmcm 

(b) 

u 

a  278  cm  s“*  =  10  km  hr 

(c) 

3 


hiqxoved  values  for  these  parameiets  wiO  be  avaibbie  when  the  measured  data  are  consolidated. 
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f. 

(d) 

»4.6x  10~^s(fore  =  8(XX)cm^s“^) 

(c) 

=  2.8  X  10~^  cm  (for  c  =  80(X)  cm^  s“^) 

(0 

« 10”^  g  cm~^ 

(8) 

,  _-2 

r^  s  10  cm 

(h) 

(3.55) 

Inserting  these  values  into  Eqs.  (3.43)  and  (3.45)  gives 

/g  =  2.6  m 

(3.56) 

/j  =  0.49  km  . 

(3.57) 

For  the  conditions  of  Eq.  (3.55),  the  small  value  of  Is  irulicates  that  turbulent  shear 

coaguladcm  will  not  be  a  factor  In  transport  since  Fc  is  essentially  constant  at  F^  1  for  all 
values  of  q  of  interest  Moreover,  it  is  hard  to  imagine  a  set  of  realistic  parameters  that 
could  change  the  foregoing  ccmclusion. 

On  the  other  hand,  Eq.  (3^7)  shows  that  /i  is  in  the  kilometer  range.  Consulting 
the  tqrptoximations  and  the  unknown  fisctors  diat  have  gone  into  the  calculation  of  h*  it  is 
not  inc(HK:eivable  that  k  could  be  as  high  as  10  km  under  some  set  of  circumstances.  This 
could  have  an  impact  cm  the  selection  and  interpretation  of  data  used  for  transport 
predictions.  Turbulent  coagulation  will  have  an  effea  on  the  prediction  and  interpretation 
of  oonoentraticm  if  Fc  varies  significantly  over  the  range  of  interest 

The  precise  criterion  about  what  constitutes  an  accq)table  degree  of  spatial  variaticm 
remains  to  be  defined  fcv  the  Kuwaiti  Oil  Fire  Field  Experiment  (KOFFE).  From  the 
analytical  fcnm  of  Fc  (Eq.  3.49),  there  will  be  some  minimum  value  of  q,  defined  as  qmin. 
beyond  which  variations  of  Fc  with  q  will  be  very  smaU.  Moreover,  since  the  dependence 
of  Fc  cm  qo  cmly  affects  the  magnitude  of  concentration  through  the  factor  exp[-//qo],  it 
follows  dutt  qmin  will  depend  only  on  /. 

Jcrfm  Gmkayne  (Ref.  17)  has  sug^ted  a  possible  criterion  for  determining  qmin. 
Cockayne's  criterion  is  that  turbulent  coagulation  will  not  be  important  in  the  range 
q^qrninifwehavetheconditicm 
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^  0.9  , 


(3.58) 


^c(*l|iiin) 


where  qmin  is  the  distance  beyond  which  /h  Fc  changes  by  less  than  10  percent  over  q^in. 
Using  Eq.  (3.49)  in  Eq.  (3.58),  we  easily  deduce 


1.1  ^  e 


+  ^<i  . 
^in 


(3.59) 


which  leads  to  die  result 

q„^S5/.  (3.60) 

Thus,  if  /  s  /i  s  0.49  km  is  used  from  Eq.  (3.57),  we  have  qmin  ==  2.45  km.  This 
of  qmin  indicates  that  turbulent  inertial  coagulation  ceases  to  be  important  for  the 
average  single  well  only  2.45  km  downstream.  On  the  other  hand,  if  /i  is  on  the  order  of 
10  km,  then  turbulent  inertial  coagulation  could  be  very  iiiqiortant  at  the  formation  of  the 
supeiplume. 

In  the  analytical  model  of  this  study,  we  have  reduced  die  dependence  of  the  spatial 
variation  of  turbulent  coagulatitm  to  the  two  parameters  qo  and  /i.  Recognizing  that  there  is 
stxne  uncertainty  in  the  numerical  values  of  these  quantities,  it  is  instructive  to  visually 
examine  the  dependence  of  rc(q)  (Ri  qo  and  h.  This  will  provide  a  preliminary  indication  of 
the  range  of  values  for  qo  and  li  that  may  be  important  in  assessment  of  transport  issues. 
The  results  are  shown  in  Figures  in-4  and  in<5.  All  of  these  curves  begin  at  q  »  qo. 

Figure  111-4  shows  the  behavior  of  Fc  as  a  function  of  q  for  assumed  values  of 
qo  -  0.1  km  and  /i  »  0.01,  0.1,  and  1.0  km,  respectively.  As  predicted,  increasing  the 
value  of  li  for  a  fixed  value  of  qo  leads  to  increasing  reduction  of  the  asymptotic  value  of 
Fc.  When  (Ij/qo)  =  10,  the  asymptotic  value  of  Fc  is  less  than  10^  of  its  initial  value.  The 
case  for  qo  =  0.1  and  /i »  10.0  is  not  shown  because  the  excessive  reduction  in  aerosol 
density  for  this  set  of  values  diminishes  interest  in  this  calculation. 

Hgure  in-5  shows  the  behavior  of  Fc  as  a  function  of  q  for  assumed  values  of 
qo  s  1.0  km  and  /i  =  0.1,  1.0,  and  10.0  km.  These  results  scale  similarly  to  those  in 
Hgure  in-4  since  the  ratios  of  Ij/qo  are  the  same  in  each  case.  In  addition,  the  respective 
asyi^ptotic  values  of  Fc  -  exp  -  [/lAiol  are  also  equal.  However,  as  shown  in  the  latter  set, 
the  scale  is  stretched  oitt  and  leads  to  signifk»nt  changes  in  density  extending  into  the  tens- 
of-ldlometers  range. 
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Figura  IIM.  Fe  M  a  Function  of  q  for  Qq  ■  0.1  km  With  it  at  a  Paramatar 


Figure  111-5.  Te  aa  a  FuncUon  of  q  tor  qo  >  141  km  With  ij  as  a  Paramatar 


IV.  CONCLUSION 


In  this  paper,  we  focus  on  the  KOF  of  1991  to  investigate  one  of  several  possible 
mass-changing  processes  that  can  limit  the  applicability  of  the  current  generation  of 
transport  codes,  i.e.,  turbulent  coagulation.  The  purpose  of  this  paper  is  to  make  a 
preliminary  assessment  of  the  conditions  where  turbulent  coagulation  is  or  is  not  an 
important  factor.  Qeaily,  if  aiibome  particles  continue  to  coagulate  in  appreciable  amounts 
during  the  transport  phase,  they  cannot  be  unequivocally  "tagged,”  and  this  potentially 
diminishes  die  validity  of  a  transport  calculation. 

The  approach  is  to  develop  an  analytical  formulation  of  the  spatial  behavior  of 
aerosol  density  as  the  particles  are  carried  along  the  turbulent  wind  field  while 
simultaneously  undergoing  coagulation.  Using  approximations  for  the  shape  of  the  aerosol 
radial  size  distribution  based  on  the  large  Canadian  forest  fire  of  1950,  we  have  been  able 
to  compute  the  time  and  space  behavior  of  aerosol  density  fOT  turbulent  shear  coagulation 
and  turbulent  inertial  coagulation.  The  latter  is  nearly  two  orders  of  magnimde  more 
effective  than  the  former  for  the  conditions  of  die  KOF. 

For  a  single  source,  the  analysis  shows  that  the  upstream  density  can  be  written  as 
N  =  N1N2,  where  Ni  is  that  part  due  to  ordinary  atmospheric  transport  processes  and  N2  is 
that  part  attributed  to  coagulation.  N2  depends  on  the  source  of  emission,  windspeed, 
atmospheric  turbulence,  and  the  characteristics  of  the  coagulation  process  and  behaves  as 
exp(--/i/q),  where  li  is  die  characteristic  dimension  associated  with  coagulation  and  q  is  the 
upstream  distance. 

Using  approximate  numbers  for  the  conditions  of  the  KOF,  /j  can  range  up  to 
several  kilometers  and,  in  some  cases,  may  be  as  large  as  10  km.  However,  in  the 
computation,  there  are  a  large  number  of  unknown  factors  that  can  be  determined  from  the 
experimental  data.  More  accurate  data  would  improve  the  estimate  of  the  coagulation 
distance  and,  hence,  define  more  precisely  the  regions  where  a  self-consistent  set  of 
measurements  can  be  integrated  with  valid  transport  models. 

The  theoretical  models  for  the  turbulent  collection  kernels  are  based  on  spherical 
shapes  for  the  aerosol  particles.  In  oil  foes,  however,  the  soot  particles  are  usually  not 
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^helical  but  are  more  strand-like.  Perh2q>s  these  particles  can  be  better  described  using  the 
theory  of  fractals. 

Preliminary  evaluaticm  suggests  that  strand-like  particles  tend  to  interact  more 
aggressively.  In  turn,  this  interaction  would  lead  to  a  larger  collection  kernel  that  would 
manifest  itself  in  a  lai^e  /j  and  would  make  turbulent  coagulation  an  important  factor  in 
defining  the  spatial  regions  where  transport  codes  can  be  c(»npaied  with  theoretical  models. 

Improvements  to  this  initial  theory  should  include  fractal  considerations  and 
improved  approximations  for  the  aerosol  size  distribution  function,  turbulence  levels, 
inhomogeneides  widiin  the  plume,  and  collecdon  efficiency. 
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APPENDIX  A 

CALCULATION  OF  MOMENTS  OF  NORMALIZED 
PARTICLE  RADIUS  FUNCTION  WHEN  s  =  3 

When  r^j » the  normalizing  constant.  A,  is  given  by  the  approximation 


For  s  s:  3,  we  obtain 


and,  hence. 


A  =  (s-l)ij;* 


A  =  2if 


g(r)  =  . 


(A.1) 


(A.2) 


(A.3) 


Using  the  fact  that  r^  »  r^  gives  the  following  results  for  the  first  four  moments; 


*u 

f  =  A  J  ^  =  2r 

?  =  a/  f  = 

h.  ^ 


(A.4) 


(A.5) 


h, 

'u 

/-A  /  rdr  =  ^lj  . 


(A.6) 


(A.7) 


The  values  of  tl  =  2  x  10-^  cm  and  nj  =  2  x  10-2  cm  are  used  in  all  of  the 
calculations. 
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